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[57] ABSTRACT 

Apparatus for characterizing the magnetic field of a 
device under test includes a device for nullifying the 
ambient magnetic field in a test environment area with 
a constant applied magnetic field, a device for rotating 
the device under test in the test environment area at 
least nearly through a complete circle, a device for 
sensing the magnetic field at a sensor location which is 
along the circumference of rotation as a function of the 
rotation angle for successive ones of a set of predeter- 
mined orientations of the device under test to obtain a 
profile of the magnetic field at the sensor location as a 
function of the rotation angle for each one of the prede- 
termined orientations, a memory for storing the pro- 
files, and a processor coupled to the memory for char- 
acterizing the magnetic field of the device from the 
magnetic field profiles thus obtained. 

18 Claims, 5 Drawing Sheets 
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PRECISION MEASUREMENT OF MAGNETIC 

CHARACTERISTICS OF AN ARTICLE WITH 

NULLIFICATION OF EXTERNAL MAGNETIC 

FIELDS 5 

ORIGIN OF THE INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract, and is sub- 
ject to the provisions of Public Law 96-517 (35 USC 
202) in which the Contractor has elected not to retain 
title. 

BACKGROUND OF THE INVENTION 

15 

1. Technical Field 

The invention is related to the fabrication of vehicles 
such as spacecraft having on board extremely sensitive 
magnetometer experiments for sensing distant magnetic 
fields whose measurements are susceptible of being 20 
distorted by local magnetic fields of components com- 
prised within the spacecraft. More specifically, the in- 
vention is concerned with ways of fabricating such 
spacecraft in such a way as to at least nearly eliminate 
such local magnetic fields of the individual components 25 
on board the craft. 

2. Background Art 

Measuring distant magnetic fields in an exploratory 
vehicle such as a spacecraft or satellite requires a sensi- 
tive magnetometer on board the craft. The problem is 30 
that a magnetometer sufficiently sensitive to measure 
distant magnetic fields (of other planets comets or aster- 
oids, for example) is sensitive to local magnetic fields as 
well. Various components aboard the spacecraft (such 
as communication equipment, on-board experimental 35 
apparatus, guidance and navigation hardware and the 
like) can have current-generated magnetic moments, 
serving as local magnetic field sources. Such local field 
sources distort the measurements of distant magnetic 
fields by the magnetometer. 

One solution has been to estimate the various local 
fields surrounding the magnetometer and compensate 
for them in processing the measurement data generated 
by the magnetometer. However, such a solution is 45 
fraught with uncertainties, not the least of which is the 
uncertainty involved in measuring or estimating the 
local magnetic field produced by a given component 
aboard the spacecraft. It would be much more desire- 
able to measure distant magnetic fields with the magne- 5Q 
tometer in the absence of local magnetic fields, thereby 
avoiding the necessity of having to “back-out” their 
distortions in the magnetometer data. However, this has 
not seemed possible because the components required 
to control a spacecraft necessarily have their own mag- 55 
netic moments. 

Accordingly, it is an object of the invention to pro- 
vide a way of quickly and accurately measuring the 
magnetic field of each individual component to be 
placed in the spacecraft. 60 

It is a further object of the invention to provide a way 
of nullifying the magnetic field of each component so 
that, when placed on the spacecraft, it will not disturb 
the on-board magnetometer. 

It is another object of the invention to provide a way 65 
of quickly measuring the magnetic field in three dimen- 
sions of each component to be placed on a spacecraft, 
permanently nullifying the component’s magnetic field 
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and then re-measuring the component’s magnetic field 
to confirm that it has been effectively nullified. 

It is a still further object of the invention to assemble 
a spacecraft almost entirely of individual components 
whose magnetic fields have been thus measured and 
permanently nullified. 

DISCLOSURE OF THE INVENTION 

The invention is a method for characterizing the 
magnetic field of a device under test, whose steps are 
nullifying the earth’s magnetic field in a test environ- 
ment area, placihg the device under test in the test envi- 
ronment area, rotating the device under test in the test 
environment area through a complete circle, sensing 
and recording the magnetic field at a point along the 
circumference of rotation of the device as a function of 
the angle of rotation of the device, re-orienting the 
device to one of a set of predetermined orientations and 
repeating the rotating and sensing steps, so as to obtain 
a profile of the magnetic field for each one of the prede- 
termined orientations, and characterizing the magnetic 
field of said device from the magnetic field profiles 
obtained for the set of orientations. In one embodiment, 
the set of orientations are the three orthogonal orienta- 
tions of said device. In another embodiment, the set of 
orientations are a set of m rotations of the device about 
an axis normal to the axis about which the rotating step 
rotates the device, where by to provide magnetic field 
profiles of said device along m great circles. In this 
latter embodiment, the characterizing step solves La- 
place’s equation for a magnetic potential along said m 
great circles simultaneously. 

The invention is also a method for fabricating and 
assembling magnetic free components, in which the 
foregoing steps are followed by the following steps for 
nullifying the magnetic field of the device under test: 
determining from the characterizing step the field 
strength and particular orientation of a permanent mag- 
net to be attached to the device so as to nullify the 
magnetic field of said device, and attaching a permanent 
magnet having the appropriate field strength to the 
device in accordance with the particular orientation, 
whereby to produce a device having virtually no mag- 
netic field at a particular distance. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is described in detail below with refer- 
ence to the accompanying drawings, of which: 

FIG. 1 is a simplified diagram of apparatus embody- 
ing the invention; 

FIG. 2 is a diagram of a series of positions with re- 
spect to a magnetic probe of a device under test whose 
magnetic field is to be characterized by the apparatus of 
FIG. 1; 

FIG. 3 is a graph of all of the successive positions to 
which the device under test in FIG. 2 must be moved in 
one embodiment of the invention; 

FIG. 4 is a flow diagram of the overall operation of 
the computer in the apparatus of FIG. 1; 

FIG. 5 is a flow diagram of one aspect of the com- 
puter operation of FIG. 4 for nulling the earth’s mag- 
netic field in the apparatus of FIG. 1 prior to measuring 
the magnetic field of a device under test; 

FIG. 6 is a flow diagram of another aspect of the 
computer operation of FIG. 1 for performing magnetic 
field mapping of the device under test; and 

FIG. 7 is a flow diagram of a further aspect of the 
computer operation of FIG. 4 for performing near field 
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measurements of the device under test using the appara- 
tus of FIG. 1. 

MODES FOR CARRYING OUT THE 
INVENTION 

Referring to FIG. 1, a turn table 10 rotates a device 
under test 12 with respect to a single axis magnetic 
probe 14. The device under test is typically a compo- 
nent of a spacecraft. The axle 16 of the turn table 10 and 
the probe 14 are stabilized on a stationary platform 18. 
A computer 20 controls a driver motor 22 connected to 
the turn table axle 16 through interface electronics 24. 
The computer 20 senses through the interface electron- 
ics a turn table position signal generated by an optical 
shaft encoder 26 on the turn table axle 16. The com- 
puter 20 also senses through the interface electronics 24 
the output of the magnetic probe 14. As the turn table 
10 rotates the device under test 12 through 360 degrees, 
the computer 20 stores in its memory 20tf a profile of the 
output of the magnetic probe as a function of angular 
position of the device (as sensed by the shaft encoder 
26). 

In order to map the three-dimensional magnetic field 
of the device under test 12, computer 20 obtains a pro- 
file of the magnetic probe output for each one of three 
orthogonal orientations of the device 12, for a total of 
three profiles. Thus, one profile of the magnetic probe 
output is obtained by the computer 20 over a 360 degree 
rotation of the table 10 with the face 31 of the device 12 
facing “up”. Then, the device 12 is re-oriented by 90 
degrees and the same step is repeated with the face 32 of 
the device 12 facing “up”. Finally, the device 12 is again 
re-oriented and the same step is repeated with the face 
33 of the device 12 facing “up”. 

Having stored the three profiles of the magnetic 
probe output in its memory 20 a, the computer 20 then 
determines from the stored data the magnetic field di- 
rection and strength of the device 12. 

The next step is to nullify the magnetic field of the 
device 12. For this purpose, a permanent magnet 36 is 
selected from a collection 38 of permanent magnets of 
different field strengths or the magnet can be manufac- 
tured in real-time to the specified field strength. In ac- 
cordance with the invention, the selected magnet 36 has 
a field strength equal to that of the device 12. The se- 
lected magnet 36 is mounted with epoxy on the surface 
of the device 12 at an orientation such that its field 
points in a direction exactly opposite to that of the 
magnetic field of the device (as determined by the com- 
puter 20). The field of the permanent magnet is thus 
equal and opposite to the field of the device, nullifying 
the device field. In this manner, a majority of the com- 
ponents to be placed aboard a spacecraft can be manu- 
factured so as to be free of magnetic fields. An entire 
spacecraft can be manufactured of components which 
have had their intrinsic magnetic fields nullified in the 
foregoing manner. 

In order to avoid mischaracterizing the magnetic 
Field of the device under test 12 due to the presence of 
the earth’s magnetic field (or any other stationary mag- 
netic field source), the environment of the device under 
test 12 must be free of any external magnetic fields. For 
this purpose, Helmholtz coils 40, 42, 44 are place along 
three orthogonal axes (labeled X, Y and Z in FIG., 
centered at the device 12. The Helmholtz coils 40, 42, 
44 consist of pairs of windings 40 a, 40 b; 42a, 42b: 44 a, 
44 b lying along a respective one of the X, Y and Z axes 
and disposed symmetrically with respect to the geomet- 
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ric center of the device under test 12. The computer 20 
sets and maintains the current in each one of the wind- 
ings so as to precisely nullify the earth’s magnetic field 
(and any other stationary magnetic field sources) in the 
5 region of the device under test 12. 

Preferably, the computer 20 is connected to a key- 
board 46 to provide manual operator control, to a moni- 
tor 48 for enabling the operator to follow the actions of 
the computer 20 and to a plotter 50 and printer 52 
10 which generate graphs and documents memorializing 
the results of the magnetic field measurements per- 
formed by the apparatus of FIG. 1. In the preferred 
embodiment, the computer 20 provides the operator 
with a pop-up menu windows environment on the mon- 
15 itor 48 which enables the operator to easily select vari- 
ous operations to be performed by the apparatus of 
FIG. 1, such as nullifying the earth’s magnetic field or 
performing magnetic mapping of the device under test 
12 or determining the magnetic field of the device under 
20 test by analyzing the three-dimensional profile of the 
output of the magnetic probe 14 obtained during mag- 
netic mapping. 

In the preferred embodiment of the invention, a more 
accurate near-field measurement of the device under 
25 test 12 is performed by obtaining the profile of the mag- 
netic probe output for a 360 degree rotation of the de- 
vice under test 12 for each one of 6 (or more) acute 
orientations of the device under test, as illustrated in 
FIG. 2. In the embodiment of FIG. 2, each acute orien- 
30 tation of the device under test 12 is achieved by rotating 
the device 12 by 30 degrees about the Y axis to the new 
orientation. Upon each 30 degree rotation of the device 
12 about the Y axis, the device 12 is held stationary on 
the turn table 10 and the turn table 10 is then rotated 
35 about the Z axis 360 degrees while the computer 20 
records the profile of the magnetic probe output as a 
function of the angular position about the Z axis. The 
foregoing is repeated for each one of six orientations of 
the device 12 displaced by 30 degrees about the Y axis. 
40 The six orientations are labelled in FIG. 3 as A, B, C, D, 
E and F. Thus, the computer 20 records a total of six 
profiles of the magnetic probe output. As will be de- 
scribed below, the computer 20 computes from the six 
magnetic probe profiles the magnetic dipole moment of 
45 the device under test with greater precision than that 
achieved in the magnetic mapping process described 
hereinabove with reference to FIG. 1. 

FIG. 4 illustrates the overall operation of the com- 
puter 20. In step 54 of FIG. 4, the computer 20 charac- 
50 terizes the magnetic probe 14 so that it can translate the 
output signal from the interface electronics 24 stimu- 
lated by the magnetic probe 14 into a magnetic field 
strength in appropriate units (such as nano-Teslas, for 
example). In the preferred embodiment, the computer 
55 memory 20 a stores a set of look-up tables characterizing 
each of many common magnetic probes by part number 
or name. In the probe characterization step 54, the com- 
puter displays a request to have the operator enter the 
probe part number or name (or both) on the keyboard 
60 46, and the computer 20 then determines the appropri- 
ate scaling factor and units to be applied to the output of 
the magnetic probe 14 received through the interface 
electronics 24. 

The next step in the process of FIG. 4 is to null the 
65 earth’s magnetic field (step 56 of FIG. 4). In this step, 
the device under test 12 is removed from the turn table 
10 and a magnetic probe is placed at the center or origin 
of the X-Y-Z axes. In one embodiment, this is accom- 
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plished by simply moving the magnetic probe 14 to the 
center of the turn table 12. The probe 14 is a single axis 
magnetic probe. First, the probe 14 is aligned parallel 
with the X axis and the computer 20 samples the probe 
output and computes the magnetic field strength along 
the X axis. From the computed X axis magnetic field 
strength, the computer 20 determines in a straightfor- 
ward manner the direction and magnitude of the cur- 
rent in the X axis Helmholtz coils 40 a, 40 b required to 
nullify the X axis measured field. The probe 14, still 
located at the origin of the X-Y-Z axes or center of the 
turn table 10, is then aligned parallel with the Y axis and 
the foregoing process is repeated to determine the re- 
quired current flow in the Y axis Helmholtz coils 42 a, 
Mb to nullify the Y axis magnetic field. The same pro- 
cess is repeated to determine the required current in the 
Z axis Helmholtz coils 44a, 44 b. The computer 20 then 
causes the interface electronics 24 to generate and main- 
tain the required current flow in each of the Helmholtz 
coils in accordance with the foregoing determination. 

The next step (step 58 of FIG. 4) is for the computer 
to notify the operator via the monitor 48 that the earth 
magnetic field nullification process is complete so that 
the apparatus is ready for a device under test to be 
placed at the center of the turn table 10. The computer 
20 offers a menu on the display 48 which enables the 
operator to choose between the magnetic mapping pro- 
cess or the near field measurement process (step 60 of 
FIG. 4). 

The computer then performs the magnetic mapping 
process (step 62) or the near field measurement process 
(step 64) depending upon the choice entered by the 
operator in step 60, Each of the mapping or measure- 
ment steps 62, 64 is followed by a unique analysis step 
(steps 66 and 68 respectively) in which the computer 20 
determines from the magnetic probe profile data the 
magnetic field of the device under test. The analysis 
steps 66 and 68 are each performed in a different man- 
ner. The result of either one of the analysis steps 66 and 
68 is that the computer 20 provides the operator with 
the magnitude of the magnetic field of the device under 
test along each of the X, Y and Z axes for a particular 
orientation of the device 12. 

The next step (step 70 of FIG. 4) is to null the mag- 
netic field of the device under test 12. This may be done 
in either one of two ways. In the preferred embodiment, 
the device field is nulled permanently by attaching a 
permanent magnet to the device, the permanent magnet 
being aligned in the direction of the device field and 
having a permanent field equal in magnitude to the 
device field. In another embodiment, the computer 20 
alters the current in the Helmholtz coils 40, 42, 44 so as 
to temporarily null the magnetic field of the device 12, 
for analytical purposes. 

In the preferred embodiment, the computer 20 deter- 
mines and in a straightforward manner how to orient 
the permanent magnet and where it should be attached 
to the device under test 12 for a maximum effective 
cancellation. The computer 20 displays the desired ori- 
entation of the permanent magnet to be attached to the 
device 12, the X, Y and Z coordinate of where it should 
be attached and the required field strength of the perma- 
nent magnet 36 (see FIG. 1). After the operator has 
successfully attached the appropriate permanent mag- 
net 36 to the device 12, the computer 20 can perform 
the final step. 

The final step (step 72 of FIG. 4) is to repeat the 
previous measurement and field determination steps 


6 

(either steps 62 and 66 or steps 64 and 68) to confirm 
that the device field has indeed been effectively nulled. 

The process for nulling the earth’s magnetic field 
performed by the computer 20 in step 56 of FIG. 4 is 
5 illustrated in detail in FIG. 5. This process is performed 
in the absence of the device under test 12. The first step 
(step 74 of FIG. 5) is to place a probe (such as the probe 
14 of FIG. 1) at the center of the turn table 10, which is 
the origin of the X-Y-Z axes. Then the computer 20 sets 
10 an axis index ”i” to 1, corresponding for example to the 
X axis (step 76 of FIG. 5). Then, the probe 14 is aligned 
parallel to the X axis, either by a servo (not shown in 
FIG. 1) controlled by the computer 20 or by requesting 
the operator to align it manually in a message displayed 
15 by the computer 20 on the monitor 48. Once the com- 
puter 20 is informed that the probe is aligned as desired, 
the computer 20 samples the magnetic probe output 
through the interface electronics 24 (step 78 of FIG. 5). 
The computer deduces the X axis component of the 
20 earth’s magnetic field from the sampled probe output 
and from this computes the amount of current required 
in the X axis Helmholtz coils 40 a, 40 b to null the mea- 
sured field (step 80 of FIG. 5). The computer 20 then 
causes the interface electronics 24 to apply the required 
25 current to the X axis Helmholtz coils 40 a 40 b (step 82). 
The probe 14 is again sampled (step 84) and if the X-axis 
field component amplitude is below a predetermined 
threshold (YES branch of step 86), then the index i is 
incremented by 1 (step 88), so that i = 2 corresponding 
30 to the Y axis, and the process returns to step 80 to repeat 
the foregoing steps for the Y axis. This includes re- 
aligning the probe 14 from the X axis to the Y axis. 
Otherwise (NO branch of step 86) the foregoing steps 
are repeated for the X axis until the proper estimate of 
35 the required null current for the X axis Helmholtz coil 
and the X-axis field component is reduced below the 
predetermined threshold, before proceeding to the Y 
axis. The foregoing process is repeated a second time 
for the Z axis (for which the index i3=3), requiring that 
40 the probe 14 be re-aligned from the Y axis to the Z axis. 
After the field components of all three axes have been 
measured nulled (YES branch of step 90), the process of 
nulling the earth’s field is finished. 

FIG. 6 illustrates in detail the process of magnetic 
45 mapping performed by the computer 20 in step 64 of 
FIG. 4. In preparation for performing this process, the 
probe 14 must be removed from the center of the turn 
table 10 and placed to one side of the table 10 in the 
location illustrated in FIG. 1. Then, the computer 20 
50 sets an index i= 1 (step 92 of FIG. 6). The computer 20 
then has the device under test 12 placed on the center of 
the turn table 10 in a first orthogonal orientation (i.e., 
with the device face 31 facing “up” along the Z axis). 
This is accomplished either by a servo (not shown) 
55 controlled by the computer 20 or by a requested to the 
operator displayed by the computer 20 on the monitor 
48 (step 94 of FIG. 6). The computer 20 then commands 
the driver 22 to rotate the turn table 10. For each rota- 
tion of the table 10 through a predetermined incremen- 
60 tal angle A 6 (step 96), the computer 20 samples the 
output of the magnetic probe 14 through the interface 
electronics 24 and stores the probe output amplitude in 
the memory 20 a (step 98). The computer 20 then inter- 
rogates the output of the optical shaft encoder 26 
65 through the interface electronics 24 to determine the 
total angle 0 through which the table 10 has rotated. If 
this amount is less than 360 degrees (NO branch of step 
100), then the computer returns to step 96 and keeps 
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repeating the intervening steps until the table 10 has 
rotated 360 degrees. Once the table 10 has rotated by 
360 degrees (YES branch of step 100), the computer 
increments the index i so that i = 2 (corresponding to the 
Y axis) and returns to step 94. This time, step 94 is per- 5 
formed by facing the next face 32 of the device under 
test 12 upward along the Z axis. The remaining steps are 
then performed as before. Then, the index i is incre- 
mented to i = 3 corresponding to the Z axis (step 102), 
and the same steps are repeated for the Z axis in the 10 
same manner. After the magnetic field components 
have been characterized for all three axes (YES branch 
of step 104), the magnetic mapping process of FIG. 6 is 
finished. 

FIG. 7 illustrates in detail the more accurate near 15 
field measurement process performed by the computer 
20 in step 62 of FIG. 4. The position of the probe 14 is 
as shown in FIG. 1. The computer 20 sets an azimuth 
index m= 1 (step 106 of FIG. 7). In step 108, the device 
under test 12 is first oriented on the turn table 10 at an 20 
orientation corresponding to the position A of FIG. 3, 
which corresponds to an azimuth angle <J> = 0 or azi- 
muth index m=l. The computer 20 then rotates the 
table 10 and for each incremental rotation through an 
angle A 6 (step 110) samples and records the output of 25 
the probe 14 (step 112). This process is repeated until 
the table 10 has rotated through 360 degrees (YES 
branch of step 114). Then, the azimuth index m is incre- 
mented (step 116) and, provided m is not greater than 
nw*(NO branch of step 118), the device under test 12 30 
is re-oriented (in step 108) by rotating it about the Y axis 
by an angle A0 = 180/m mfl *to the next position (position 
B of FIG. 3). The measurements of steps 110 through 
114 are then repeated. The foregoing process is re- 
peated until all values of the azimuth index m (up to and 35 
including m ma x ) for the positions C through F of FIG. 

3 have been characterized (YES branch of step 118). 

Step 68 of computing the magnetic field of the device 
under test from the magnetic mapping data produced by 
the magnetic mapping process of FIG. 6 is performed 40 
by the computer 20. The computer 20 simply deter- 
mines from each of the three profiles of the magnetic 
sensor output (corresponding to the three orthogonal 
orientations of the device) the angle at which the maxi- 
mum magnetic probe output was obtained, for a total of 45 
three angles (one for each orthogonal orientation). For 
each of the three angles, the computer 20 stores the 
magnitude of the corresponding magnetic probe output 
and determines the corresponding magnetic field 
strength. By treating these three angles' as Euler angles, 50 
the computer 20 uses straightforward spherical trigo- 
nometry to compute the field strength components 
along the X, Y and Z axes of the device 12 for a given 
device orientation (e.g., where the device face 31 faces 
“up” along the Z axis of FIG. 1. 55 

In the preferred embodiment, the computer 20 deter- 
mines from the X, Y and Z components of the measured 
device field the field strength and orientation of a per- 
manent magnet to be mounted on the device under test 
12 so as to nullify the magnetic field of the device under 60 
test. 

Step 66 of computing the device magnetic field from 
the near field measurement data produced in the near 
field measurement process of FIG. 7 employs a unique 
application of Maxwell’s equations. Specifically, the 65 
computer uses the near field test data taken for each 
value of the azimuth index m to solve the matrix equa- 
tion (9) of Appendix A hereof. The angle 4> of FIG. 1 is 
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the angle <f> of equation (3) of Appendix A. The azimuth 
index m of FIG. 7 is the index m of equation (3) of 
Appendix A. In the preferred embodiment, the azimuth 
index m assumes the values 1, 2, 3, 4, 5 and 6 so that the 
rotations of FIG. 3 are performed at 30 degree intervals, 
thus providing magnetic measurements in six great cir- 
cles. All of magnetic field strengths deduced from the 
six rotations performed in the near field measurements 
process of FIG. 7 are substituted appropriately into the 
matrix equation (9) of Appendix A by the computer 20. 
The computer 20 solves for the tensor quantities Aj 1 , 
Bi 1 , Aj° of the matrix equation (9) of Appendix A, 
which are the magnetic dipole moments of the device 
under test 12 along the X, Y and Z axes, respectively. 
From this, the computer 20 determines the radial and 
angular components H r and He of the device magnetic 
field as 

H r — 2 M r ~ 3 cos 0 and 
Hq—M sin $, 

where r is the distance from the dipole to the point at 
which it is measured, M is the magnitude (absolute 
value) of the dipole moment vector and 0 is the angle 
between the dipole moment vector and the radius vec- 
tor r between the dipole moment and the point of mea- 
surement, in accordance with well-known principles. 
From this, the computer 20 determines in straightfor- 
ward manner the field strength and orientation of the 
permanent magnet to be attached to the device under 
test to nullify its magnetic field at a predetermined dis- 
tance. 

Industrial Applicability 

The invention is useful in assembling exploratory 
vehicles such as spacecraft having a highly sensitive 
on-board magnetic sensor with magnetic-free compo- 
nents whose magnetic fields have been nulled by first 
characterizing the . magnetic field of each component 
thereof so as to not affect the on-board sensor. 

While the invention has been described in detail by 
specific reference to preferred embodiments thereof, it 
is understood that variations and modifications thereof 
may be made without departing from the true spirit and 
scope of the invention. 

Appendix A 
Copyright 1990 
Nearfield Analysis Theory 

Before the data is submitted to the analysis it is ad- 
justed for three things First, offsets which occur be- 
tween curves due to external fields or test equipment 
adjustment are removed. Since the values at the poles 
should be equal on all six traces, each curve is adjusted 
by the average difference between the north and south 
poles of all six traces. Thus the values at the poles on all 
six traces are made approximately equal. Second, the 
average value of all data points is subtracted from every 
data point. This is because the divergence of the field 
over the spherical surface is zero (Maxwell). Third, 
each data point is weighted according to the relative 
spherical surface area which it represents. Points near 
the poles are closer together and represent less area in 
the upcoming surface integration. 

The analysis itself is based on Maxwell’s Second 
Law, which states that the curl of a magnetic field on a 
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closed surface is zero. Thus the field outside the mea- 
surement sphere can be expressed as a scalar potential 
that is a solution of Laplace’s equation. 


dr 2 


10 

Cjj are the normalized trigonometric coefficients, 
obtained by substituting r, 6 , and <f> for each data 
point into equation (3) 

N is the total number of points (NC*NP) 

The error between the calculated and measured value 
is given 


2 c 2 ct> cotfl cQ 1 ? 2< t> _ q 

r or + p. jl dB r^sin 2 # d<b 2 

10 

<t>(0,6) = ( 2 > 

2 — ^ — * 2 [A n m cos md> + B n m sinm<f>] P n m (cosO) 

n=l r*+ ] m = 0 

D (3) 15 

Br = ~ ~ = 


35 (6) 

ej S= ''2 ^ ( CijAi — B r j) 

N N 35 (?) 

2 e/= 2 2 -B ri ) 

i- 1 

To minimize the amount of error to find the least 
squared best fit to equation (7). 


2 — J-T 2 [A n m cosmd) + B n m sinm(J>] P n m (cos0) 
«= 1 m— 0 


N . 

A* 


A 7 35 , 

= 2 2 {CjjAj — B ri ) 

i=\j=\ 


( 8 ) 


where, 

4> is the magnetic scalar potential 


20 If equation (8) is differentiated to find local minima and 
then expanded and expressed in matrix form: 


, (9) 

(2cy j)/Jj + (2c,‘.ic,-.2)^2 + . . . + (2c/.}C/.35M35 = 2ci A Brj 

aci.2Ci.0Ai + (2c7 2 )^2 +■••+ (2C/.2C/, 35)^35 = 2ci,\Br I 


(2c/.35C/,iM35 -+■ (2c/. 35C/.2>^2 + • • • + (2^.35^25 = 2Ci,35i?r/ 


(cc ] )A = cB r 


B r is the radial magnetic field (calculated) 
r is the radius of the measurement sphere (probe dis- 
tance) 

0 is the table angle (point) 40 

<b is the tilt angle (curve) 

P n m is the Schmidt function of degree n and order m 
(scaled associated Legendre functions) 
n is the degree of the spherical harmonic term 
m is the order of the spherical harmonic term 4 5 

A n m and B n m are the spherical harmonic coefficients 
(multipole moments) 

With the experience obtained through the years at 
JPL, it has been determined that the highest practical 
degree of moment needed is n = 5. It has also been found 5 q 
that the number of curves (NC) required to achieve an 
accurate measurement is six. The number of points per 
trace (NP) to be used in the analysis is selected by the 
operator at the time of analysis. With these quantities, 
equation (3) can be represented for each point as a sum- 55 
mation of 35 terms. 

B ri = C U \A i°+ C a A i 1 + 0,3*1 1 + C/,4^2°+Q5^2 1 - 
+ 0,6*2* +♦ . . + 0, 3si?5 5 

B ri = C it \A i + C it2 A 2 + 0,3^3 +. • • + O, 3*435 
When shown for all data points equation (4) becomes: 


If the c matrix is inverted the unknowns can be found. 

A = (cB)*(cc) 00) 

The multipole moments can be extracted from these 
unknowns as follows: 
monopole moment: 

A 0 °=0 

dipole moment: 

Al 1 =Dj: 

Bi 1 = D P 
Ai° = D z 

quadrupole moment: 

Qli = 3* A 2 2 -A 2 ° 

Q 22 =-3* A 2 2 -A 2 ° 

Q33-2 A 2 ° 

Ql2 = 3* A 2 2 
Qi3 = 3* A 2 1 
Q23=3^ B 2 i 

Coding of Important Functions 

datgen: averages and equates poles 
ggen: calculates area weighting factors 
pnm: calculates associated Legendre functions, given: n, 
m, and cos0; uses mathematical definition - differenti- 
ates cos 2 0 — 1 

spnm: scales associated Legendre functions to obtain 
Schmidt polynomials 

smfcn: calculates normalized trigonometric coefficients 
regan: calculates Schmidt coefficients from equation 
(9b), uses matrix inversion 

crintm: calculates dipole and quadrupole moments from 
Schmidt coefficients 


A 7 35 

Br = il Xj lx Ci ri 


(5) 


65 


where 

Ay are the unknowns A m n and B m n 


(4) 60 
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senfld: calculates field at alternate coordinates due to 
earlier calculated multipole moments 


12 

spacecraft project. The system only uses them mathe- 
matically when conducting nearfield analysis. 


Smoothing 


Derivation of DUT Angles 


Transient noise appears superimposed on top of the 
true measured data waveforms. This noise is a result of 
60 Hz wiring in the vicinity, movement of metal ob- 
jects. variations in the earth's field, etc. The Schonstedt 
magnetometer has two rear panel controls which ena- 
bles internal filters. The first is a power line filter to 
keep 120v AC noise from reaching the data. It should 
always be kept ON. The second filters the probe read- 
ings. and is adjustable to 10 Hz, 100 Hz, or 1 kHz. This 
knob should be kept on 10 Hz to attenuate as much 
noise as possible. 

The rest of the unwanted noise which passes through 
these filters can be attenuated with software. This soft- 
ware smoothing is best understood when considered in 
the frequency (f) domain. The noise which is filtered 
out is of a high frequency, but the true data is usually 
low frequency. Thus, a low-pass filter, chosen with the 
correct cutoff frequency (fc), will eliminate the un- 
wanted noise, leaving the desired data. 

In signal processing, the data waveform is convoluted 
with a window' function {w(f)} to perform this filtering. 
Windows with sharp edges produce unwanted “sine” 
sidelobes during convolution (Gibbs effect). The best 
window to minimize this phenomenon is the Blackman 
Window', defined as: 

*</) = 0.42 + 0.5 *cosin *f/f 0 ) + 0.08 *cos(2it *f/f Q ) 

The cutoff frequency is chosen differently, depending 
on where the smoothing is requested. If you are prepar- 
ing data for input to nearfield analysis: 


/o - 


Number of Points per Input Curve 
360 


cycles/point 


At any other time, such as during a test, or while exam- 
ining: 


f 0 — 10/360=0.2778 (cycles/point) 

Maximum Radial Zero-to-Peak Field 

Maximum radial zero-to-peak field at one meter is 
calculated during analysis of magnetic mapping data- 
files. It can only be computed if all three rotations were 
performed at the same probe distance. The algorithm 
used: 


* (B x 2 +B } ?+bM 

where, 

B r is the maximum radial zero-to-peak field computed 
at one meter (nT) 

r is the probe distance (m) 

Bjt, By, B z are the maximum peak-to-peak fields of 
each rotation (nT) 

example) A magnetic mapping is performed at 1 m. 
The pk-pk fields for all three axes are measured to 
be 25 nT. The resultant field at 1 m is calculated to 
be 15.31 nT. 

Derivation of DUT Coordinates 

The DUT coordinates specify where, in XYZ space- 
craft coordinates, the center of the DUT is located. The 
values are reported in meters and are defined by the 


5 The angles w'hich associate the fixed spacecraft coor- 
dinate system ( X s / C , Y s /c, Z s / C ) to the DUT measure- 
ment coordinate system (Xz>,Y/), Z o) are termed Z-Y-Z 
Euler angles. The procedure for obtaining these angles 
(a,/?,0 is listed below: 

10 1. Begin with the two frames coincident. 

2. Rotate the DUT frame clockwise, (as seen from the 
spacecraft coordinate origin), about the Z s / C axis by 
an angle a. 

3. Rotate the DUT frame clockwise about the Y s / C 
I 5 axis by an angle /?. 

4. Rotate the DUT frame clockwise about the Z s / C 
axis by an angle T. 

5. The DUT coordinate system is now oriented in its 
measurement position. 

20 

Units of Magnetic Field Strength 

A magnetic field surrounds a magnetic source. Its 
strength is measured in units of nanoTesla (nT). A nano- 
Tesla is also frequently called a gamma. Another com- 
25 mon unit of field strength is the Gauss (G), w'hich is also 
called the Oersted The relationships between these units 
are given in Table E-l. 


TABLE E-l 



Magnetic Field Conversion Factors 



Equivalent field in units of: 

Field 

nanoTesla 

Gauss 

Tesla 

1 nanoTesla 

1 

10“ 5 

io- 9 

(gamma) 
1 Gauss 

10* 

1 

10- 4 

(Oersted) 
1 Tesla 

10 9 

10 4 

1 


Units of Magnetic Moments 

A magnetic moment is a measure of the strength of a 
magnetic source. Its value is independent of reference 
frame. The chief unit which is used is “Gauss centimeter 
cubed” (G*cm 3 ), also called the “pole centimeter” (po- 
45 le*cm) Other units are “Ampere meter squared” 
(A*m 2 ), and “Weber meter” (Wb*m). Conversion fac- 
tors for these units are listed, in Table E-2. 

TABLE E-2 

Magnetic Moments Conversion Table 


Equivalent moment in units of: 

Moment G * cm 3 nT * m 3 A * m- Wb • m 


1 G * cm 3 

1 

0.1 

10 ~3 

1.26 X 

io- 9 

(pole * cm) 
1 nT * m 3 

10 

j 

10-2 

1.26 X 

10~ 8 

55 1 A * m 2 

10 3 

10 2 

1 

1.26 X 

10“ 6 

1 Wb * m 

7.94 X 30 8 

7.94 X I0 7 

7.94 x 10 5 

1 



Relating Magnetic Fields and Moments 

60 The relationship between a magnetic moment and the 
field which surrounds it will vary depending on the 
type of source. Simple dipoles, multipoles, long wires, 
Helmholtz coils, etc all have different relationships. 
They are too complex to describe in detail here. 

65 The relationship between a pure dipole moment and 
its magnetic field is given by: 


H r ~ 2 M r 3 cos& 
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Hq—M r~ 2 sinQ 

where, 

H r is the component of the magnetic field which is 
directed along r (Gauss) 

H 74 is the component of magnetic field which is 
perpendicular to the radial magnetic field, and in 
the plane of M (Gauss) 

r magnitude of the vector r between the dipole mo- 
ment and the point where the field is measured 
(meters) 

M magnitude of the dipole moment vector M 
(G*cm 3 ) 

0 is the angle between M and r. 

At a distance of 1 m, a dipole that is aligned with the 
DUT X-axis, with strength 500 G*cm 3 , has a peak 
radial field of 10“ 3 Gauss, or 100 nT. 

What is claimed is: 

1. A method for characterizing the magnetic field of 
a device under test, comprising: 

nullifying the ambient magnetic field in a test envi- 
ronment area, said nullifying comprising: 

(a) measuring and retrievably storing three orthog- 
onal components of the ambient magnetic field 
strength in said test environment area, 

(b) computing in processor means a requisite cur- 
rent to be provided in respective orthonal Helm- 
holtz coils for nullifying the corresponding or- 
thogonal magnetic field strength components, 
and 

(c) maintaining a constant current equal to said 
requisite current in each respective Helmholtz 
coil whereby to nullify said ambient field so as to 
refrain from changing said constant current 
while performing subsequent sensing steps: 

placing the device under test in the test environment 
area; 

rotating the device under test in the test environment 
area at least nearly through a complete circle; 

sensing and recording the magnetic field at a sensor 
location which is along the circumference of rota- 
tion of the device as a function of the angle of 
rotation of the device; 

re-orienting the device to a successive one of a set of 
predetermined orientations and repeating said ro- 
tating and sensing steps, so as to obtain a profile of 
the magnetic field at said sensor location as a func- 
tion of the angle of rotation of the device for each 
one of said set of predetermined orientations of said 
device under test; 

characterizing the magnetic field of said device from 
the magnetic field profiles obtained for said set of 
orientations. 

2. The method of claim 1 wherein said set of orienta- 
tion comprise three orthogonal orientations of said de- 
vice. 

3 . The method of claim 1 wherein said set of orienta- 
tions comprise a set of m rotations of said device about 
an axis normal to the axis abut which said rotating step 
rotates said device, where by to provide magnetic field 
profiles of said device along m great 

4 . The method of claim 3 wherein said characterizing 
step comprises solving Laplace’s equation for a mag- 
netic potential along said m great circles simulta- 
neously. 
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5 . The method of claim 1 further comprising the 
following steps for nullifying the magnetic field of said 
device under test: 

determining from said characterizing step the field 
5 strength and particular orientation of a permanent 

magnet to be attached to said device so as to nullify 
the magnetic field of said device; and 

attaching a permanent magnet having said field 
strength to said device in accordance with said 
10 particular orientation. 

6 . The method of claim 4 wherein said solving com- 
prises solving Laplace’s equation with Legendre poly- 
nomials whereby to obtain the magnetic quadrapole 
moment of said device under test. 

15 7. A method of assembling an exploratory vehicle 

having an on-board magnetic field sensor whereby indi- 
vidual components of said vehicle introduce virtually 
no magnetic fields to said on-board magnetic field sen- 
sor, said method comprising the following steps carried 
20 out for each one of said components prior to their inclu- 
sion in said vehicle: 

nullifying the ambient magnetic field in a test envi- 
ronment area, said nullifying comprising: 

(a) measuring and retrievably storing three orthog- 

25 onal components of the ambient magnetic field 

strength in said test environment area, 

(b) computing in processor means a requisite cur- 
rent to be provided in respective orthogonal 
Helmholtz coils for nullifying the corresponding 

30 orthogonal magnetic field strength components, 

and 

(c) maintaining a constant current equal to said 
requisite current in each respective Helmholtz 
coil whereby to nullify said ambient field; 

35 placing the one component in the test environment 
area; 

rotating the one component in the test environment 
area at least nearly through a complete circle; 

sensing and recording the magnetic field at a sensor 
40 location which is along the circumference of rota- 
tion of the one component as a function of the 
angle of rotation of the one component; 

re-orienting the one component to a successive one of 

- a set of predetermined orientations and repeating 
45 said rotating and sensing steps, so as to obtain a 
profile of the magnetic field at said sensor location 
as a function of the angle of rotation of the one 
component for each one of said set of predeter- 
mined orientations of said one component; 

50 characterizing the magnetic field of said one compo- 
nent from the magnetic field profiles obtained for 
said set of orientations; 

determining from said characterizing step the field 
strength and particular orientation of a permanent 
55 magnet to be attached to said one component so as 
to nullify the magnetic field of said one component; 
and 

attaching a permanent magnet having said field 
strength to said one component in accordance with 
60 said particular orientation, whereby to produce a 
component for inclusion in said vehicle which at 
least nearly does not introduce a magnetic field to 
said on-board sensor. 

8. The method of claim 7 wherein said set of orienta- 
65 tions comprise three orthogonal orientations of said one 

component. 

9. The method of claim 7 wherein said set of orienta- 
tions comprise a set of m rotations of said one compo- 
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nent about an axis normal to the axis about which said 
rotating step rotates said one component, where by to 
provide magnetic field profiles of said one component 
along m great circles. 

10. The method of claim 9 wherein said characteriz- 5 
ing step comprises solving Laplace’s equation for a 
magnetic potential along said m great circles simulta- 
neously. 

11. The method of claim 7 further comprising the 
following steps for nullifying the magnetic field of said 
one component: 

determining from said characterizing step the field 
strength and particular orientation of a permanent 
magnet to be attached to said one component so as 15 
to nullify the magnetic field of said one component; 
and 

attaching a permanent magnet having said field 
strength to said one component in accordance with 
said particular orientation. 

12. The method of claim 10 wherein said solving 
comprises solving Laplace’s equation with Legendre 
polynomials whereby to obtain the magnetic quad- 
rapole moment of said device under test. 

13. Apparatus for characterizing the magnetic field of 25 
a device under test, comprising: 

means for nullifying the ambient magnetic field in a 
test environment area, comprising: 

(a) a set of three orthogonally oriented Helmholtz 
coils placed near said test environment area, 

(b) means for measuring three orthogonal compo- 
nents corresponding to said three Helmholtz 
coils of the ambient magnetic field strength in 
said test environment area, 

(c) means responsive to measurements obtained 
from said means for measuring for computing a 
requisite current to be provided in respective 
orthogonal Helmholtz coils for nullifying the 
corresponding orthogonal magnetic field 
strength components, and 

(d) control means responsive to said means for 
computing current for maintaining a constant 
current equal to said requisite current in each 
respective Helmholtz coil whereby to nullify 
said ambient field so as to refrain from changing 
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said constant current during subsequent sensing 
steps; 

means for rotating the device under test in the test 
environment area at least nearly through a com- 
plete circle; 

means for sensing the magnetic field at a sensor loca- 
tion which is along the circumference of rotation of 
the device as a function of the angle of rotation to 
the device for successive ones of a set of predeter- 
mined orientations of said device with respect to 
said means for sensing so as to obtain a profile of 
the magnetic field at said sensor location as a func- 
tion of the angle of rotation of the device for each 
one of said set of predetermined orientations of said 
device under test; 

means for storing the profiles obtained by said means 
for sensing; and 

processor means coupled to said means for storing for 
characterizing the magnetic field of said device 
20 from the magnetic field profiles obtained for said 
set of orientations. 

14. The apparatus of claim 13 wherein said set of 
orientations comprise three orthogonal orientations of 
said device. 

15. The apparatus of claim 13 wherein said set of 
orientations comprise a set of m rotations of said device 
about an axis normal to the axis about which said rotat- 
ing step rotates said device, where by to provide mag- 
netic field profiles of said device along m great circles. 

16. The apparatus of claim 15 wherein said processor 
means comprises means for solving Laplace’s equation 
for a magnetic potential along said m great circles si- 
multaneously. 

17. The apparatus of claim 13 further comprising: 
35 processor means coupled to said memory for determin- 
ing from the output of said means for characterizing the 
field strength, a location and an orientation of a perma- 
nent magnet to be attached to said device so as to nullify 
the magnetic field of said device and means for display- 

40 ing said field strength, location and orientation of said 
permanent magnet to be attached. 

18. The apparatus of claim 16 wherein said means for 
solving comprises solving Laplace’s equation with Le- 
gendre polynomials whereby to obtain the magnetic 

45 quadrapole moment of said device under test. 

* * * * * 
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